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I. INTRODUCTION 

There have been suggestions that mass, production and decay properties of some of the isoscalar scalar mesons are 
consistent with a situation where a glueball should have mixing with quarkonia states in a mass range roughly below 



2 GeV. The idea that a glueball configuration is present in the scalar meson spectrum was suggested in [l| and later 
extensively used as for example in Refs. The discussion on possible evidence for the scalar ground state glueball 

has dominantly centered on the scalar-isoscalar resonances /o(1500) and /o(1710). The claim that at least partial 
glueball nature can be attributed to these states is to some part based on the mass predictions of Lattice QCD. In 
the quenched approximation the lightest glueball state is predicted to be a scalar with a mass of about 1650 MeV 
^ and uncertainty of around 100 MeV. But note that first results in an unquenched calculation could indicate a strong 
downward mass shift towards 1 GeV [8[ . For reviews on the experimental situation of scalar mesons and their possible 
y—i . structure interpretation see for example @. 

ly-^ ' An important object for studying the nature of the scalar resonances above 1 GeV has been the strong decay 
patterns observed for these states. Thereby several analyses conclude that the three scalar states /o(1370), /o(1500) 
and /o(1710) result from a mixing of the glueball with the nn — \/l/2 (uu + dd) and ss states of the scalar 3 Pq 
quarkonium nonet. Furthermore, the nonappearance or the strength of the appearance of these scalar states in 77 or 
in central pp collisions have been argued to give further signals for a possible quantification of the glueball component 
residing in the scalars @-@|. A study of the strong and radiative three-body decays of J/tp with J ftp —> V + PP(yy) 
and J/tp — > 7 + PP(VV), where V and P are vector and pseudoscalar mesons, opens a further opportunity to 
test the nature of the scalar mesons / (1370), /o(1500) and / (1710). The strong three-body J/tp decays with 
PP = ir+ir- and K + K~ in the hnal state have been studied by the Mark III [10j, DM2 [ll[ and BES II [H, [H| 
collaborations. Particularly the BES II [13j experiment with a much larger statistics indicated a clear signal for 
/o(1370) — > in the cf>ir + ir~ data. An enhancement in the cpKK data can be fitted by interference between the 

/o(1500) and the /o(1710). In [h34l6| BES II also reported on partial wave analyses on radiative decays of the type 
J/tp — > 7PP and J/tp — > yVV. In particular, they indicated the branching ratios B(J/tp — > 7/0 (1710) —> yKK) |14| . 
B(J/ip -> 7/o(1710) 77T+7T-), B(J/ip 7/o(1500) ->• 77T+7T-) [l| and B{J/ip -> 7 / (1710) yuico) [Ig. 

Theoretical studies of J/tp decay into a vector meson and two pseudoscalars have already been done in Refs. [l7T - fl9l ] . 
The impact of the effects of coupled-channel dynamics has been investigated in pj} and has then been improved on 
the basis of chiral perturbation theory (ChPT) and unitarity constraints in fl8j . A coherent study of the production 
of the mesons /* = /o(1370), / (1500) and / (1710) in J — > Vft -> VPP {V = <p,w) has been considered in [| in the 
framework based on a glueball-quarkonia mixing scheme involving the scalar mesons fi. An analysis of the role of 
the scalar mesons /o(1370) and /o(1710) in J/ip radiative decays has been performed in Ref. [20T [2l|. 
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It was shown in Refs. [22|,[23[ that the dominant contribution to the Okubo-Zweig-Iizuka (OZI) suppressed radiative 
decays of the J/tp comes from the photon emission from the initial state. This means that the strong J/tp decays 
involving isoscalar vector mesons (like u> and (p) are not related to the radiative decays of the J/tp via the vector 
meson dominance (VMD) mechanism. For this reason the branching ratios of strong and radiative J/tp decays are 
compatible. The three-body J/tp decays can be factorized [3j in terms of the two-body rates. Therefore, the analysis 
of the J/tp three-body decays can shed light on the two-body transitions of the J/tp and fi states: J/ip — > fcV, 
J/tp — > /,7, fi — > V"f and /; — > 77. The electromagnetic decays of the scalar mesons fi — > V7 and fi — )• 77 have been 
studied in detail using chiral approaches, nonrelativistic and light-front quark models (see e.g. Refs. @, Hj-[28j|). 

The purpose of the present paper is to analyze the role of the scalar mesons fi in the strong and radiative three- 
body decays of J / tp using a chiral Lagrangian approach suggested and developed in [f| H^, [30] • In these works we 
originally studied the strong and electromagnetic decay properties of scalar mesons above 1 GeV. The isoscalar scalars 
were treated as mixed states of glueball and quarkonia configurations. Later on we extended the phenomenology to 
the study of decay properties of excited tensor, vector and pseudoscalar mesons. Present evaluation is meant as a 
continuation of the full decay analysis presented in jBj but now we include the additional constraints set by J/tp 
decays involving these scalars. At this stage we do not consider explicitly the mixing of the scalar resonances induced 
by the continuum decay channels, we also neglect final state interaction in the decay modes (note that inclusion of 
final-state interaction in J/tp three-body decays has been done in [lj|) _ both interaction mechanisms are higher-order 
effects in the framework of our perturbative considerations. We are interested to give predictions for scalar resonance 
contributions to J/tp decay modes where the physical nature of the scalar resonances fi — Jo (1370), Jo (1500) and 
/o(1710) can possibly be tested. 

In the present paper we proceed as follows. In Sec. II, we present the effective Lagrangian which will be used 
for the calculation of the resonance contributions of scalar mesons fi to the matrix elements of strong and radiative 
three-body decays of the J/tp. In Sec. Ill we present our results for the resonance contributions of the scalar mesons 
fi to the decay widths of the J/tp. A short summary is given in Sec. IV. 

II. APPROACH 

In Refs. @, H^, [3(| we presented the lowest-order chiral Lagrangian describing strong and radiative decays of 
pseudoscalar, scalar, vector and tensor mesons. This Lagrangian was motivated by chiral perturbation theory [3ll — 
l34j . Here we extend the formalism by including the additional couplings J/tpfiV, J/tpfij, fiVV and fiV^j. These 
interaction terms are necessary for elaborating the contributions of the scalar mesons fi to the strong and radiative 
three-body J/tp decays. The full Lagrangian relevant for the J/tp meson decays involves the J/tp meson, the scalar 
glueball G, the nonets of pseudoscalar V, scalar S and vector V mesons with 

£ cff = ^(D ll UDi>rf + X +) + \(D P SD»S - MlS 2 ) + ^Gd^G - M 2 G G 2 ) 

- \& »W» v V»W pv - hl 2 W^Wn - \{d^ v dPJ pv - -Mjj^JV") 

c 9 , „, ... c! 



+ c s d (Su^) + cl n {S X+ ) + ^C(u/) + ^=G(X+) + C v mix + Ct 
-9 c 



+ <f e (SF+F+» v ) + ^=G(F+F+n + ~J-(S{W»\F+}) + -±G(W^F+) (1) 

+ c^SW^Wn + ^=G(W^Wn + c s h J^(SW^) + c s w J^{S){W^) 

+ ^J^GiW^) + J^F^ (c n k N + 4S + 4G) . 

The symbols (..) and { } occurring in Eq. ([l]) denote the trace over flavor matrices and the anticommutator, re- 
spectively. \N) = I (uu + dd) /y/2) , \S) — \ss), \G) are the nonstrange, strange quarkonia states and the glueball, 
respectively @. The constants c] describe the couplings of scalar fields and of the glueball to pseudoscalar, vector 
mesons, J/tp and to photons. For the vector mesons we use the tensorial representation in terms of antisymmetric 
tensor fields W^ v (SU(3) nonet of vector mesons) and J pv (J/tp meson) [see details in [34], l3q]. The couplings of 
the scalars to vector mesons and photons {c s el c 9 e , Cy, (A, c s v and c^) can be constrained by vector meson dominance 



3 



(VMD) and SU(3) flavor symmetry as: 



„2 „2 



-■' - ^cg. (2) 



V6 

As we already stressed in the introduction, the couplings of J/ip to scalars and photons contained in c\ (i = n,s,g) 
are not constrained by VMD and, therefore, are independent on the couplings c s h and c 9 h . In the derivation of the 
expressions in Eq. we also used SU(3) relations for the vector-meson to photon couplings: 

■9p7 = 3guy = ~^a<t>i (3) 

with g pl = 0.2 fixed from data on T(p° e + e~) [36]. The other notations are standard for the basic blocks of the 
ChPT Lagrangian [3l]-[34j]: U — u 2 — exp(iPv / 2/ F) is the chiral field collecting pseudoscalar fields in the exponential 
parametrization, F is the pseudoscalar meson decay constant, and V denote the chiral and gauge-invariant 
derivatives acting on the chiral fields and other mesons. Furthermore, x± = u*x u ± u\^u with x = 2B(s + tp), s = 
M + ... and F/f v = u'F^vQu + uF^Qw, where F pv is the stress tensor of the electromagnetic field. The charge 
and the current quark mass matrices are represented by Q = ediag(2/3, —1/3, —1/3) and M = diag(™, m, m s ) where 
m = (m u + md)/2 and m s are the nonstrange and strange quark masses. In the following we restrict to the isospin 
symmetry limit m u = ra^. The terms >Cj^i X arid £f lix give rise to flavor singlet-octet mixing in the pseudoscalar 
and scalar sector, while for the scalar case quarkonia-glueball mixing is included in addition [see details in Q]. In 
particular, the mass matrix involving nonstrange \N) and strange \S) = \ss) quarkonia states and the glueball \G) 
has the form [|| 

/ M 2 N fV2 e \ 
ML C = fV2 Ml f , (4) 
\e f Ml) 

where / and e are free parameters controlling the quarkonia-glueball mixing scheme. The orthogonal physical states 
are obtained by diagonalization of M^ are with the transformation matrix B as 

BM 2 m . c B T = M) = diag(M / 2 i , M% , M% ) . (5) 

The eigenvalues of Mj represent the masses of the physical states /i = /o(1370) , f% = /o(1500) and fs = /o(1710) . 
These diagonal states with i = fi, / 2 , /3, are then given in terms of the bare states as 

(6) 

The J/tjj three-body decays which proceed through the scalar resonances are described by the diagrams shown 
in Figs. 1-4. The diagram in Fig.l contributes to the strong decay J/ip — > V fi — > VPP. There are two graphs 
in Fig. 2 relevant for the radiative decay J/ip — > Vjj: one described by the two-step process J/ip — > fiV — > 
[Fig. 2(a)] and the other one with J/ip — > fjpf — > ^77 [Fig. 2(b)]. Fig. 3 contains the graph responsible for the radiative 
decay J/ip — > fci — > ^PP involving two pseudoscalars in the final state. Finally, the diagrams of Figs. 4(a) and 
4(b), representing the two-step transitions J/ip — > ff) — > jVV and J/ip /jV — > "fVV, contribute to the decay 
J/ip — > jVV with two vector mesons in the decay products. The corresponding matrix elements and decay widths 
resulting from these specific graphs are given in full detail in the Appendix. Note, that the three-body decay widths 
of J/ip are factorized in terms of two-body decays of the subprocesses J/ip — > fiV, J/ip — > fij, fi —> PP, fi — > VV, 
fi -t V-f and fi -> 77. 

III. RESULTS 

A. Mixing schemes of the scalar mesons /o(1370), ,fo(1500) and /o(1710) 

In Ref. [6] we already presented a detailed analysis of strong and radiative decays of the scalar mesons fi = 
/o(1370), /o(1500) and /o(1710) e.g. considering different scenarios of quarkonia-glueball mixing. Bare masses, 
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mixing parameters and decay couplings of the quarkonia and glueball configurations were extracted from a fit to 
experimental masses and well-established decay rates into two pseudoscalar mesons. In the present discussion we also 
indicate a general fit of these parameters, especially for the decay constants c d and of the direct glueball decay. 
In Ref. Q we constrained these parameters by either studying the flavor-symmetry limit - c 9 n = or by matching 
these constants to first lattice results for the glueball decay. In the following fits no constraints are set on these 
decay parameters, but previously deduced mixing schemes will approximately result again from this procedure. The 
quantities entering in the fitting procedure and the results obtained from both fits are displayed in Table I. In the fit 
here we also include the additional decay mode /o(1500) — > r\rj . As in [y] (T/ 3 )2P is the sum of partial decay widths 
into two pseudoscalar mesons approximated by the total width. By releasing the constraints on the direct glueball 
decay we obtain two solutions where the quality of fits in terms of xtot ^ s somewhat improved compared to the ones 
of Ref. @- We consider two options resulting from local minima in Xtot'- Scenario I, the mass of the bare glueball G is 
on the lower end of the unquenched lattice predictions with tuq ~ 1.5 GeV; Scenario II, the mass of the bare glueball 
G is larger with mc ~ 1.7 GeV. The fit parameters for these two schemes are fixed as follows: 
Scenario I: 

M N = 1.485 GeV , M G = 1.482 GeV , M s = 1.698 GeV 

/ = 0.068 GeV 2 , e = 0.236 GeV 2 , c s d = 8.8 MeV (7) 
c s m = 2.2 MeV , c 9 d = 1.8 MeV , c 9 m = 27.7 MeV 

Scenario II: 

M N = 1.360 GeV , M G = 1.686 GeV , M s = 1.439 GeV 

/ = 0.23 GeV 2 , e = 0.30 GeV 2 , c s d = 6.5 MeV (8) 
c s m = 5.5 MeV , c 9 d = - 2.0 MeV , c 9 m = 48.3 MeV 

The corresponding mixing matrices B\ (Scenario I) and B\\ (Scenario II) are: 

/ 0.75 0.6 

Bi = -0.59 0.8 -0.14 | , (9) 
V-0.29 -0.05 




/0.79 0.29 0.55 \ 
Bu = 0.57 - -0.82 . (10) 
\0.23 -0.96 0.17 / 

For solution I the glueball component dominantly resides in the /o(1500) (with Xtot — 20), while for the second 
one (with xtot — 15) /o(1710) is identified with the glueball. Both level schemes of the bare states and mixing 
scenarios contained in Bjjj are very similar to the original solutions of Ref. [6j, where the direct glueball decay 
is included (denoted in @ as third and fourth solution). Releasing the constraint on the direct glueball decay 
parameters allows a further fine-tuning of the results. The two solutions also represent the situation discussed in the 
literature, when studying the mixing of the glueball with quarkonia in the /o(1370), /o(1500) and /o(1710) sector 
- the glueball resides dominantly in the /o(1500) or the /o(1710). Qualitatively both solutions are in acceptable 
agreement with the experimental data, except for the underestimate of ^k^(i430)-^Ktti which as in the previous work 
Q remains unresolved. In the context of the effective chiral theory the second solution is slightly favored. We present 
additional theoretical ratios not contained in Table I, compared for example to experimental data from the WA102 
Collaboration pzj (see Table II). Note that in contrast to the result of T }l ^ RK /Y = 0.46 ± 0.19 of the WA102 
Collaboration 37] the analysis of the OBELIX Collaboration [H] results in 0.91 ± 0.20. 



B. Hadronic and radiative J/i/> decays 

In this subsection we analyze the role of the scalar resonances fi in the hadronic J ftp decays J/ip — > Vfi — > VPP. 
Details of the calculations are presented in the Appendix: analytical formulas for the matrix elements and the 
decay rates. Our strategy for determining the parameters involved in these processes is the following: using data 
(central values) on the three-body decay modes J/ip — > (4>, w)/o(1710) — > (cp,uj)KK [3(| (see Table III) as well as the 
predictions of our approach for the two-body decay modes of the scalars (see Table I) we fix the couplings c s h , c s w , c g h 
in the effective chiral Lagrangian ([1]) as 

= 1.092 x 10~ 3 GeV" 1 , c s w l = -0.626 x 10~ 3 GeV" 1 , cf = 10~ 4 GeV" 1 , 
cf 1 = 1.340 x 10~ 3 GeV" 1 , c s ^ u = -0.924 x 10~ 3 GeV" 1 , cf; 11 = 0.722 x 10~ 3 GeV" 1 . (11) 



5 



The indices I, II refer to the respective mixing scenarios. With these couplings we determine the strong couplings 
9JfiV [they are linear combinations of c| and c 9 h (|A6[) ] and the electromagnetic couplings gj/ i7 of the physical /j 
states. In Table IV we list the results for the effective couplings gjhv and <7j/ i7 and the corresponding two-body 
decay widths. 

In Table III we give our final results for the hadronic three-body decays of J/tp for both scenarios (I and II). 
Results are compared to data (3(| and the predictions of Ref . Q , where a similar mixing scheme [l| has been used to 
determine B. Especially in scenario I we obtain a reasonable description of present data on the hadronic J/tp decays, 
while in version II some discrepancies occur. In particular, one can see that in scenario I we can reproduce the central 
value of the ratio 

= Br( J/tP -> M -> cpKK) _ 4 

Br(J/V> f 3 0J ojKK) 3 1 ' 

while in scenario II this quantity is close to 2. The explanation for this quantitative difference between the two mixing 
schemes is quite simple. The ratio R is approximately equal to the ratio of the corresponding strong two-body decay 
modes of J/ip, which is expressed in terms of the elements of the mixing matrix B as: 



Br(J/V -> huj) 



zl(B 31 /V2-B 



:s:S I 



c s h B 33 + c 9 h B 32 /V3 + c* w (B 31 V2 + B 33 ) 



(13) 



The ratio depends crucialy on the value of the combination of mixing matrix elements with A = B 3 \/\2 — B 33 . In 
scenario I we have A = —1.16 and therefore it is possible to generate the observed value for the ratio R with the 
appropriate choice of free parameters cf, c? and c s w (as we did in our fitting procedure). In case of scenario II this 
quantity has the value A = —0.007 ~ 0, and now the ratio R is always close to 2 independent on the choice of the 
relevant parameters c|, and c s w . After the mxing matrix B is constrained by a fit to the properties of the scalars 
(as compiled in Table I) scenario II is less preferable when analysing the hadronic J/tp decays. 

Analysis of the radiative J/tp decays J/tp — > 7 + PP(VV) and J/tp — » P + 77 can produce further important 
information about the electromagnetic structure of the scalar resonances involved. The relevant matrix elements 
involve further unknown couplings (cj, Cj, c s e , c 9 , <£, cf, c£, c|, c 9 k ) from the effective Lagrangian ([T]). Note, only 
5 parameters from this set are independent, because the parameters c s j, c s e and c?, cf are expressed through the 
parameters c s v and c 9 via VMD constraints ©. 

The other parameters involved were already fixed from the previous analysis of the strong J/tp decays. We find 
that acceptable results for the radiative J/tp decays can be achieved with the following choice of the parameters , 

r a r " r s „9 . 
c vt c k> c fe) c k- 



c s / = 2.815 GeV" 1 , c 9 ' 1 = 0.138 GeV -1 , 

c/f = 0.241 x 10~ 2 GcV" 1 , cf = -0.313 x 10~ 2 GcV -1 , c 9 / = -0.271 x 10~ 3 GeV" 1 . 
cf 1 = 3.036 GeV" 1 , c 9 / 1 = 1.898 GeV" 1 , 

c n k Jl = 0.132 x 10" 2 GeV" 1 . cf 1 = 0.187 x 10~ 2 GeV" 1 , cf 1 = -0.352 x 10~ 2 GeV" 1 . 



(14) 



Knowledge of the couplings c s v and c 9 gives access to the effective couplings of scalar mesons to vector mesons and 
photons. They are related by the VMD relations ^ and ([3]). With the use of Eq. (TT4")) and the VMD relations we 
determine the effective couplings of the scalar /j mesons to photons and vector mesons as listed in Table V. 

Our predictions for the radiative decays of J j tp and the scalar mesons are given in Tables VI and VII. Note that 
we reproduce all known data for radiative decays of the J/tp involving the /j states [36| : 

Br( J/tp -> / 27 ) = (1.01 ± 0.32) x 10~ 4 , 

Br(J/^ -> / 37 7 7T7r) = (4.0 ± 1.0) x 10" 4 , (15) 
Bt(J/1> -)■ / a7 ->■ jKK) = (8.5± -J) x 10~ 4 , 
Br( J/tp -> / 3 7 -> jujuj) = (3.1 ± 1.0) x 10" 4 . 

In case of the scalar mesons we set our results in comparison to the ones of other approaches and with available 
data. In particular, the results of Ref. 24 1 and [25| are given in the form (L, M, H), where L=Light, M=Medium and 
i/=Hcavy correspond to the three possibilities for the bare glueball mass: lighter than the bare nn mass, between the 
nn and ss masses, heavier than the ss mass. The results of Ref. (2(| are listed in the form (K, P), where K = KK 
and P = 7T7T indicate the contributions of the intermediate KK and nn loop, respectively. 
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The calculated branchings of the two- and three-body decays of J/ip and the /, mesons satisfy the following 
approximate ratios: 

Br( J/ip -> fa -> ttttt) D Br(/i -> tttt) 



Br( J/V -+ fa -+ KKi) ~ Bi(fi^KK)' 

Ra Br( J/i, -+ fi V) " - Br( J/iP /i7 ) " Bf(/ * PP) " (16) 

We want to illustrate these approximate constraints in case of the scenario I. The set of the radiative three-body J/ip 
branching ratios Ru = 1.04, R21 = 3.80, R31 = 0.56 is similar to the corresponding set of the strong two-body ft 
branching ratios R12 — 0.94, R22 = 3.61, R32 = 0.43. Also, the ratios RJz PP and Rf P are approximately equal to the 
corresponding branchings of the strong two-body decays of the scalar mesons fa. E.g. in case of f 3 = /o(1710) and 
K K in the final state we have 

R§g K = 0.45 ~ R^ K = 0.44 ~ Rf 4 K = 0.45 ~ Bi(f 3 -tKK) = 0.49. (17) 

Similar consistency is found in case of the scenario II with: 

J?u=0.51 ~ i?i2=0.44, # 21 = 0.41 ~ # 21 = 0.39, R 31 = 0.64 ~ i? 31 =0.49, (18) 
R% 3 KK ~ R^ K ~Rl\ K = 0.38 ~ Br(h^KK) = 0.42. (19) 



IV. CONCLUSIONS 



In conclusion, we present an analysis of the role of the scalar resonances /o(1370), /o(1500) and /o(1710) in the 
strong and radiative J/ip three-body decays. This analysis can shed light on the nature of the scalar resonances 
fi and the possibility that these states arise from glueball-quarkonia mixing. We tested two cases for the glueball- 
quarkonia mixing schemes: scenario I, the bare glueball dominantly resides in the /o(1500); scenario II, the scalar 
/o(1710) contains the largest glueball component. We found that the first scenario is more consistent with present 
data, especially for the case of the strong three-body decays of J/ip. The results are also consistent with available 
data on radiative decays. The detailed set of predictions for the radiative decays can serve to further distinguish 
between the different structure assumptions of the scalars. It would therefore be extremely useful to have more data 
on the radiative processes. 
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Appendix A: Matrix elements and decay widths of the J/ip three-body decays 

For the calculation of the three-body decays of the J / ip we introduce several notations p is the momentum of the 
J/ip and q, qi, q2, are the momenta of the final state particles: V(q), P(qi), P{q 2 ) [for the decay J/ip — > V fi — > VPP], 
V(q), 7(91), 7(*) [decay J/iP -+ /,V( 7 ) -> Vyy], l(q),P(qi),P(q 2 ) [decay J/i\> -> f a -> 1 PP] and 7 (g), V{q x ), V(q 2 ) 
[decay J/ip -+ fiV(j) ^VV]. 

We also define the invariant variables Si(i = 1, 2, 3): 

P = q + qi + <72 , 

si = (q + qi) 2 = (p-92) 2 , 

52 = (qi + q2? = (p - qf , 

53 = (q + q2) 2 = (p-qi) 2 . (Al) 

The matrix elements describing the J/ip three-body strong and radiative decays proceeding through a scalar resonance 
/, are denoted by M^X U ly , with super- and subscripts referring to the graphs of Figs. 1-4 in an obvious notation. 
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The explicit expressions for these amplitudes read as 



Mi = M(J/> -> fiV -> V^PP) = ^(p)c^(g) ff/,pp (g^pq-p^) 



Mtt = 



1 



Mi. - s 2 - iM f T u 



(A2) 
(A3) 



M| 4 - s 2 - iM fi T fi ' 
M{J/1> -> /(V -> ^77) = e$(p) £^(9) < a (?0 ^(ft) flJ/,7 9/<Vi 

(g^pqi -p a qip) (g^p mi -qm») 



(9mpPQ2 ~ Pp<hii)(9va qqi - q a qiv)Yri 



Mj. -S 3 - iMyXy, 
1 



My 3 . - si - iM/Tj 

M ln = M(J/ip -> 7/, ; — > 7PP) = ej(p)e* l/ (q) gj/,v gnpp (g^pq-p^) 

Miv = M( J/^ /<V(7) ->■ = M#+M$, 

M^ = M(J/^^/ i 7^7^) = ^(p)< I/ (g)er(9i)e^(g 2 )5J/ <7 P/ i yv 

x (g f ,upq-Puq fl ,)(g a /3qiq2- qi02 a )jj2 nrr~F > 

My — S 2 — *Myd y i 

M$ = M(J/ib -+ UV jVV) = e%p) e™{q) e^(gi) e^(g 2 ) 



M? - s 2 - iM h T h 



(A4) 
(A5) 



(ff^P* - p/jq^Xsw. qqi - q a qi 



1 



M|. -Sa-iMyTy, 

1 



Ji,f2 



M/ 4 - S X - iM fi T fi 

where Mf i and Ty 4 are the mass and width of the scalar meson /j. Here, the coefficients gjf t v, <?/;PP) gfiW, 
9 J fill 9fiVi an< i ff/;77 are the couplings involving mixed scalars, which are related to the couplings of the effective 
Lagrangian ([lj involving unmixed states as: 



gjf^ 

9 film 

9fiKK 

9f iPP 



2B. a c s h + 2^B l2 c 9 h + 2V2{B tl V2 + B a )c s w , g JM = -2B i3 c s h - zJ^B^ - 2{B a V2 + B i3 )C(A6) 



2 Bn 
F 2 y/2 



{Ml - 2Ml)c s d + 2Mlc% 



gjfi<i> 

2 Bj2 

F 2 -s/3 



(My 2 t - 2Mj) c * + 2Mjc^), 
2 B i2 



—=c s v B lY + —=c 9 v B i2 , gy,, 



%/2 



v/3 



4 



4cf,B i3 + —c 9 B l2 , 



16/5 



g J/i7 = 2{B ll c k l + B i2 c 9 k + B l3 c s k ) , g hll = — [— =B a + B l3 c; 



32 



9 V\/2 *7 e 3V3 

/2 2 2 

-Bacy , g fi<tn = -B i3 c) + ^-j=B i2 c 9 f 



(A7) 

(A8) 
(A9) 
(A10) 
(All) 



In above expressions the Bij are the elements of the matrix B relating the mixed (/1, f 2 , f 3 ) and unmixed (N, G, S) 
scalar states as in Eq. ([6]). If the process is kinematically allowed, the couplings 3,7/^, gjft-y, 9fiPP> 9fiW, 9fiV-y and 
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5/i77 define the corresponding two-body decay widths of J/ip and the physical /j states as: 



2 



(A12) 



Vo-kM u Y M|. 



H/, ^ />/>, V,,,. I -^f , VI:-!) 



r(j/</^/ i7 ) = ^^^(^l-^J , (Ais) 

r(/i->V 7 ) = f^M^l-^) 3 , (A16) 
r(/,^77) = J« 2 4 77 M 3 . (A17) 



where A(x, y, z) — x 2 + y 2 + z 2 — 2xy — 2xz — 2yz is the Kallen function. The factor 

'§, PP = irir 
2, PP = KK 



[ 3 pp — ww 

N PP ={ 2 ' - (A18) 



takes into account the sum over charged modes: 

r(/i -> 7T7T) = T(fi -> 7T+7T-) + r(/, -> 7r°7r°) , 

r(/i->#jr) = r(/ J ^x+x-) + r(/ J ^x°/?°). (Ai9) 

Note, the couplings g/iWi 9fiVj, Qfi-y-y an d gjfoV) 9Jfi-y are constrained by VMD and SU(3) flavor symmetry relations 
(see also Eqs. © and ©): 

ff/.77 = ffp7 (^/;P7 + gff/^7 + -g-S/*^^ = 3p 7 ^/iW + g9/iWW + g9f,4><pj > ( A2 °) 

1 \/2 

5/.P7 = 9p-f9fr PP , 5/,^7 = gflWS/f^, ff/,07 = -^-9pi9h<t><t> , (A21) 

5J/ l7 = ffP7 Q^j/^ + ^-gjfi^j (A22) 

with g p7 = 0.2. 

Next we write down the expressions for the J/ip three-body decay widths. The decay width T\ = T(J/ip — > foV — > 
VPP) is given by: 

{Mj-Mvf 4,1 2 2 

where the sum is performed over the polarizations and where 

(Mj-Mvf 



\i/ 2 (M 2 ,M 2 lS2 ) AM% f x(M 2 M 2 , ]+ a M 2 M 2\ 

/. / ^ 2(M 2 _ S2)2+r2 M 2 W 1 - — [KM j, M v ,s 2 ) + 6MjMy j , 



IMt 



t = M 2 P + UM 2 + M 2 - S2 ±^ 2 { S2l MlM 2 )Jl-^\. (A24) 



S2 

Using Eqs. (|A12I) and (IA13|) we express Ti through the two-body decay widths T(J/ip -» /jV) and -> PP): 

T - r< ' J ^ ~* ^ ~* PF " > M ^ R (A25) 

Tl tt AV2(Mj, M\ , M* ) ^ - 4M| ^P^ji ^ - O + 6M 2 M 2 R ' (25) 
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The decay width Tu = T(J '/i/j — > fiV(j) — > V-f-f) is calculated according to the expression: 

(Mj-M v ) 2 

J ds 2 J ds^WnW (A26) 



r T 



1536tt 3 M 3 



pol 



where 



± Af 2 + Af 2 -s 2 X^(s 2 ,M 2 ,M v ) 
s iji — 2 2 ' 

The sum over the polarizations is rewritten as 

^2 \ M n\ 2 = 47r 2 a 2 ( g 2 JfiV g 2 fnj R x + g 2 Jhl g\ Vl R 2 + g JfiV 9jf iy 5/,y 7 5/, 77 #3 J (A28) 

pol ^ ' 

where 



/>>J = (M 2 -S2? + M*r* ( A(M "' M ^ S2) + 6m2jM v ) ' ! A2< ' 1 



R 



2 



(Mj- Sl ) 2 (M 2 - Sl ) 2 , (Af 2 -s 3 ) 2 (M 2 -s 3 ) 2 



(M 2 - si) 2 + M I Tj z (Mj. - s 3 ) 2 + M 2 T 2 
(Af 2 - si)(M 2 - s 3 ) + Mj.Tj. 

((Ml - Sl y + m 2 ^((M 2 - s 3 ) 2 + m 2 r 2 ; 



(s x s 3 - MjMy) + siMjMy , (A30) 



B - gl / (M 2 -s l)( M 2 -s 2 ) + M 2 T 2 

^ 3 " (Ml - s 2 f + M)Y\ \ (Ml — si) 2 + Af 2 .r 2 . 
(M? — S3) (Ml - S2) + M? r 2 f „ o o\ 

+ W-4' + «?n (s ' + M3<) ) ' <A31) 

The decay width r m = T(J/ip -> £7 -> 7PP) is given by: 



M 2 , s+in 



Fni 768tt 3 M 3 



where 



jj J ' 384^ 2 Af 3 ^7 V s 2 (Af 2 - s 2 ) 2 + M 2 r 2 

4M| po1 AMI J ' U U 



sf,u = M 2 P + (l ±] jl-^j. (A33) 

Using Eqs. (IA13|) and (|A15I) we again express Tni through the two-body decay widths T(J/tp —> fij) and r(/j — > PP) 
as: 

Ml 



AMI 



T{J/^^fa)T{f i ^PP)Ml f 4M 2 (Af 2 -s 2 ) 3 fA . 



7T (Afj - Af 2 .) 3 ^/M 2 - 4MJ, 



p 4M| 



s 2 (M 2 - s 2 ) 2 + m 2 r 2 . 



The decay width T IV = T{J/tj) — ?> fiV{^) — » 7UF) is determined according to: 

Mj s tiv 
J - pol 



4 AT, 2 , 
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where 



, Jlv = Mj + ^(l ± ^/l-M). (A36) 

Again, the sum can be expressed as 

|Mv| 2 = ira(gj fn gj. vv Qi + g 2 JfiV g 2 . Vl Q 2 + gjf.v g.jf n 9uw guv 1 Qs J (A37) 
pol ^ ' 

with 

= 7772 1 , , f2 p2 f4-2^(Mj + 2M2) +s 2 (M 4 + 8M 2 M 2 +6M 4 } 



(M 2 . — S2) 2 + M'fT'j, 

As 2 M 2 J M^{M' 2 j + 3M y ) + 6MjM v ) , (A38) 



° 2 = (m] - Sl y + m 2 r 2 ~ 2s ' (M ' + 2M ^° + S ' (M ' + 8M2/M ^ + 6M ^ } 

- 2siM^(Mj + 5MjM y + 2M y ) + M y (Mj + AM 2 My + My) 

+ (M 2 -s,Y+M 2 T 2 (*» ~ 2S ' (M ' + 2M v ] + S 3( M / + 8M 3 M v + 6M V ) 

- 2s 3 M y (M j + 5MjM y + 2M V ) + M y (Mj + 4MjM y + My) 



(M£ - Sl ) (Ml s 3 ) + M 2 Y\ , 2 2 



+ ((Ml- S S + M 2 r^((Ml- S3 r + M 2 r 2 J I ^ + W + ^(Mj + My) 



- 2(*i + s 3 )M y (Mj + 2My) + 2My(2Mj + 5M 2 )J , (A39) 

O 1 f ( M l- S ^ M l- S ^ + M l T l (** , o,2 Af 2, f 2 , q 2 M 2 (M 2 , ^ 

Q3 - (M 2 . - s 2 ) 2 + M 2 Y 2 , { (M 2 . — si) 2 + M 2 Y 2 ^ ^i^ + ^M^ + ^M^M. + M,) 

- 2(sj + s 2 )M 2 jMlr(M 2 j + 2My) + MjM y (Mj + 4M 2 M y + 2M y ) 



- 2(s 3 + s 2 )MjM^(Mj + 2M y ) + MjM y (M] + 4MjM y + 2M y ) . (A40) 
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FIG. 1: Contribution of the scalar resonances /; = / (1370), / (1500), / (1710) to the decay J/ip -> frV ->■ VPP 




7 V 



(a) (b) 
FIG. 2: Contribution of the scalar resonances /; = / (1370), / (1500), /o(1710) to the decay J/ip -> fiV(-y) -» V77 
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FIG. 3: Contribution of the scalar resonances /; = / (1370), / (1500), /o(1710) to the decay J/ip ->■ / l7 -> 7PP 




(a) (b) 
FIG. 4: Contribution of the scalar resonances /< = /o(1370), /o(1500), /o(1710) to the decay J/ip -> fiVfr) -»■ 7VF 



TABLE I: Masses and decay properties of the scalar mesons fi. 



Quantity 


Data [36] 


Fit I 


2 

X 


Fit II 


2 

X 


Mf 1 (MeV) 


1350 


± 


200 


1432 


0.16 


1231 


0.36 


M fa (MeV) 


1505 


± 


6 


1510 


0.72 


1510 


0.70 


Mf 3 (MeV) 


1720 


± 


6 


1720 


~ 


1720 


~ 


T f2 ^ w (MeV) 


2.07 


± 


0.87 


1.2 


0.97 


1.02 


1.45 


Tfr*m ( MeV ) 


5.56 


± 


0.98 


2.8 


7.98 


4.70 


0.76 


r h ^RK (MeV) 


9.37 


± 


1.09 


10.4 


0.96 


8.5 


0.19 


T j2 ^^ (MeV) 


38.04 


± 


2.51 


37.7 


0.02 


38.38 


0.02 


Ffa^vv/F f 3 ^KK 


0.48 


± 


0.15 


0.25 


2.45 


0.25 


2.45 




0.41 


± 


0.14 


0.43 


0.03 


0.49 


0.36 


T h 2P (MeV) 


137 


± 


8 


137 


~ 


136 


0.01 




0.35 


± 


0.16 


0.29 


0.15 


0.29 


0.13 




0.88 


± 


0.23 


0.8 


0.11 


1.06 


0.64 


rjf.(i430)^^ (MeV) 


251 


± 


74.40 


64.7 


6.27 


41.64 


7.92 


Xtot 


19.82 




14.99 



TABLE II: Comparison of strong scalar decays with WA102 data [37| . 



Quantity 


WA102 data [37] 


I 


II 


r f ± -> k k 1 r s\ ->tvtv 


0.46 ± 0.19 


1.07 


2.27 


T/i — «7>7 /T/l -urir 


0.16 ± 0.07 


0.22 


0.4 


T fl ^ 2 p (MeV) 


small 


314 


131 


^ is -» v v ' 1 ^ fa -> t t 


< 0.18 


0.89 


1.4 



TABLE III: Branchings of the J/t/> hadronic decays in units 10 



Meson 


4>KK 


Data [36] 


ojKK 


Data [36] 


07T7T 


Data [361 


CJ7T7T 


Data [361 




I 


II 


[3] 




I 


II 


[3] 




I 


II 


® 




I 


II 


® 




/3 


3.6 


2.5 


3.6 


3.6 ±0.4 


4.8 


5.2 


4.8 


4.8 ± 1.1 


1.7 


1.3 


0.40 




2.2 


2.7 


0.53 




h 


0.5 


2.3 


0.19 


0.8 ±0.5 


0.2 


1.1 


1.38 




2.0 


9.4 


0.77 


1.7 ±0.8 


0.6 


4.6 


5.60 




h 


4.1 


4.7 





0.3 ±0.3 


1.3 


8.7 


0.08 




4.3 


2.4 


0.01 


4.3 ± 1.1 


1.4 


4.4 


0.83 





TABLE IV: The branchings of the J/ip two-body decays in 10 4 and effective couplings in 10 s x GeV 



Quantity 


I 


II 


Data [36] 






- 0.60 


- 1.91 








0.57 


2.2 




9Jf3" 




- 1.60 


- 1.80 




9Jfi4> 




1.02 


1.35 








-1.00 


2.18 




9 J fat 




-1.40 


1.27 




9.1 in 




2.34 


2.12 




9Jf2~t 




- 1.53 


- 1.53 




9Jf 3 i 




- 7.42 


7.98 




Br(J/> - 


/iw) 


2.2 


27.5 




Br(J/V> - 


/2W) 


1.8 


13.6 




Br(J/V> - 


/aw) 


10.7 


13.5 




Br(J/V> - 




6.8 


14.9 




Br(J/V> - 


/a0) 


5.9 


28.0 




Br(J/^ - 




8.1 


6.6 




Br(J/V> - 


/i7) 


2.6 


2.6 




Br(J/i/) - 


M) 


1.0 


1.0 


1.01 ± 0.32 


Br(J/V - 


/s7) 


17.7 


20.5 





TABLE V: Effective couplings of ft in units GeV" 1 . 



Coupling 


I 


II 


9 fin 


0.30 


0.41 


9 fill 


- 0.21 


0.13 


9 fall 


- 0.01 


- 0.08 


9fipi ~ 3ff/i"7 


1.24 


1.60 




0.29 


0.75 


9fipi = 3(?/ 2 oj7 


- 0.90 


0.98 


9f2<t>l 


- 0.13 


- 0.94 


9f3Pl = 3<?/ 3 u>7 


- 0.47 


- 0.44 


9fs4>i 


0.10 


- 0.20 


9f\pp — J/i"" 


6.19 


8.01 


9fiU 


3.09 


8.00 


9fipp — 9feutt> 


- 4.44 


4.91 




- 1.33 


- 9.97 


9fspp ~ 3/3"^ 


- 2.35 


- 2.19 


9f 3 H> 


10.74 


- 2.16 



15 



TABLE VI: Branchings of the J/0 radiative decays in units 10 



Mode 


I 


II 


Data [36] 


J/1>- 


• hi 


7Tf7T 


1.65 




0.42 








* J27 ->■ 


7Tf7r 


0.34 




0.34 






J/4*- 


■ hi 


7Tf7T 


4.48 




5.02 




4.0 ± 1.0 


J/i>- 


' hi 


-yKK 


1.58 




0.82 






J/1>- 


- hi 


-yKK 


0.09 




0.08 






J/i>- 


■ /s7 


-yKK 


8.00 




7.85 




8-5-0.9 


J/i>- 


' /i7 




0.33 




0.22 








> /a7 -> 


JLOLU 


0.13 




0.16 






J/1>- 


■ hi 


JUILJ 


3.09 




3.07 




3.1 ± 1.0 




' /i7 ->■ 


Pll 


2.0 xlO" 


2 


1.7 xlO" 


2 




J/4>- 


' /i^(7) 


—¥ CJ77 


2.2 xlO" 


3 


3.3 xlO" 


3 




J/1>- 


• /i0(7) 


— > 077 


0.8 xlO" 


3 


2.5 xlO" 


3 






' /27 -> 


Pll 


0.9 XlO" 


2 


1.1 xlO" 


2 




J/4>- 


' /sw(7) 


— > wyy 


1.1 xlO" 


3 


1.5 xlO" 


3 




J/i>- 


' /2 0(7) 


— > 077 


0.4 xlO" 


4 


0.5 xlO" 


2 




JH>- 


' /s7 -> 


Pll 


6.0 XlO" 


2 


6.5 xlO" 


2 




j/1>- 




— > CJ77 


0.7 xlO" 


2 


0.7 xlO" 


2 




j/i>- 


• h<Pil) 


— > </>77 


0.16 




0.8 xlO" 


2 





TABLE VII: Decay widths of /, — ¥ Vj and fi —¥ 77 transitions (in keV) in comparison with other theoretical predictions. 



Mode 


Ref. [24] 


Ref. [25] 


Ref. [26] 


Ref. [28] 


Ref. [6] 


Our (I, II) in keV 


h - 


-> 77 




(1.6,3.91^,5.611;*) 




0.35 


1.31 


(11.1, 13.4) 


h- 


-> 77 








0.35 




(6.5, 2.4) 


h - 


■* 77 




(0.92,1.31^,3.011*) 




0.019 


0.05 


(0.02, 1.2) 


h - 


-> p7 


(443, 1121, 1540) 


(150,3901^,53011?°) 


(79 ± 40, 125 ± 80) 




726 


(1441, 957) 


h - 


■> W7 






(7 ±3, 128 ±80) 




0.04 


(156, 102) 


h - 


-> 07 


(8,9,32) 


(0.98, 0.83±8;1S, 4.5±tg) 


(11 ±6,-) 




0.01 


(27, 30) 


h- 


-> P7 












(990, 1209) 


h- 


■> W7 












(108, 131) 


h- 


-> 07 












(8, 447) 


h- 


-> p7 


(42, 94, 705) 


(24,55ll 6 4 ,410l^) 


(100 ±40, -) 




24 


(519, 450) 


h- 


■> W7 






(3.3 ±1.2,-) 




82 


(57, 49) 


h- 


-> 07 


(800, 718, 78) 


(450,400l^,36lH) 


(15 ±5,-) 




94 


(1300, 53) 



